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TECHNICAL MEMORANDUM X-64841 


SSME/SIDE LOADS ANALYSIS FOR FLIGHT CONFIGURATION 

INTRODUCTION 


This memorandum presents the results of a comprehensive Monte Carlo 
dynamic loads analysis of the Space Shuttle Main Engine ( SSME) for the side 
load transients produced by flow separation in the engine bell during ground 
ignition. The side load excitation is a transient phenomenon which occurs dur- 
ing ignition of an engine that has a high expansion ratio nozzle like the SSME. 
The previous publication of this document is revised herein to include a more 
refined structural model of the engine and the orbiter thrust structure and to 
include an active thrust vector control system (in the previous analysis the 
actuators were simulated with "equivalent?’ springs). The results presented 
herein pertain only to the flight configuration; test configurations that include 
constraints on the engine different from the flight configuration will require 
additional analysis. 


STRUCTURAL MODEL 


The structural modeling and vibration analysis were accomplished by the 
Structural Dynamics Section, S&E-ASTN-ADS, MSFC. The dynamic response 
of the engines is greatly influenced by the stiffness characteristics of the orbiter 
thrust structure at the interfaces with the engine. The orbiter thrust structure 
flexibility matrix was supplied by Rockwell International, Space Division. The 
flexibility coefficients are different for each engine and for each engine/ thrust 
structure interface point. This necessitates a structural model for each engine. 
The engine locations, their orientations, and coordinate axes are shown in 
Figure 1. The description of the structural models and their vibration analyses 
are documented in MSFC memorandums S&E-ASTN-ADS (74-23) and S&E- 
ASTN-ADS (74-25). The frequency of the pendulum mode (first mode) was 
determined to be 8.49 Hz, 9.42 Hz, and 8. 15 Hz for engine numbers 1, 2, and 
3, respectively. 

Figure 2 presents the NASTRAN node points and coordinate system for 
the structural idealization. Table 1 shows the degrees of freedom and grid 
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point coordinates. The structural models of all the engines are identical; only 
the backup structure differs for each engine. Two coordinate systems are 
identified in Table 1. Coordinate system no. 4 is fixed with its origin at the 
gimbal point. Coordinate system no. 3 is rotated with respect to coordinate 
system no. 4 to obtain the various start positions specified for each engine. 

Five propellant ducts were included in the structural model and are identified in 
Table 2. A sketch indicating the most pertinent nodes for this analysis is shown 
in Figure 3. It is noted that the actuators are not included as part of the struc- 
tural model. For this analysis the mode shapes of the "dangling'’ engine (engine 
without actuators) are used. The actuators are simulated as feedback control 
systems via the transient response equations (see Appendix for details). 

The structural elements monitored for loads are identified in Table 3. 
The NASTRAN sign convention is used for element loads and is shown on Figure 
4. A model damping value of Z, - 0. 01 was used in all the "dangling" modes. 
Modes were retained through 100 Hz in the transient response analysis. 

TVC ACTUATOR MODEL 


All actuators in the orbiter thrust vector control system are not identical. 
However, for this study all were assumed to be identical to the pitch actuator on 
engine 1. The actuator locations are shown in Figure 1. The pitch actuator on 
engine 1 is being developed by MSFC and presently is the only one for which a 
dynamic model is available. All the actuators are expected to behave in approxi- 
mately the same way, dynamically. A description of the actuator model is 
provided in MSFC memorandum, S&E-ASTR-GM- 140-74, SSME TVC Actuator 
Block Diagram for Transient Load Analysis, dated April 16, 1974. The block 
diagram is shown on Figure 5 and the actuator parameters are tabulated in 
Table 4. The procedure used to couple the actuator model to the structural 
model is explained in the Appendix. 


SIDE LOAD FORCING FUNCTIONS 


The SSME side load forcing functions were obtained by scaling J-2S side 
load data to the SSME. The J-2S data were gathered and analyzed by Rocketdyne 
from J-2S firing data and J-2S dynamic shake test data. Rocketdyne supplied 
an analog tape containing 24 side load transients and a 40 Hz bandwidth white 
noise signal from which the transients were obtained. These transients were 
derived by shaping the white noise signal as indicated in the following block 
diagram. 
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WHITE NOISE (0 -► 40 Hr) 




The J-2S side load data were extracted at the bell exit; i.e., grid point 33 on 

Figure 2. It is observed that the magnitude of the side load transient is directly 

related to F* , the rms force per 1. 6 Hz bandwidth of the white noise. Data 
max 

scatter was observed in determining F* from various measurements ( load 

max 

cells and accelerometers). Thus this parameter is treated as a random 
variable in the Monte Carlo analysis. The statistical distribution for F* was 
generated from test data supplied by Rocketdyne. 1 The data were taken from 
test 043 for the time slice corresponding to the operating point for maximum 
side loads. The data points included measurements from both load cells and 
accelerometers. The distribution is presented in Figure 6. Rocketdyne’s 
assessment of all valid test data had previously yielded a mean +3o- value of 
1600 lb rms per 1. 6 Hz bandwidth for J-2S side load. Correspondingly, the 
distribution in Figure 6 was adjusted to yield the same mean +3cr value. 

The analog data tape supplied by Rocketdyne was digitized by R-COMP- 
RRV at MSFC. Plots of the digitized transients are presented in Figures 7 
through 30. The digitized transients are referenced to F* max = 1432 lb rms per 

1. 6 Hz resolution bandwidth. It was observed that the digitization at 400 samples 
per second tended to generally decrease the peak value of the transient by about 
10 percent. 

Information with regard to directionality of the side load excitation was 
not available from J-2S data. For this analysis the direction of the side load 
excitation, 9 , ( Fig. l) is assumed variable from 0 deg to 360 deg and is 
considered to be constant for each side load transient. Thus, 6 is a random 
variable for the Monte Carlo analysis with uniform distribution between 0 deg 
and 360 deg. 


1. Internal Letter SSME 73-3199, Statistical Evaluation of J-2S Side Load Data 
for SSME, November 15, 1973. 
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To obtain side load excitation applicable for the SSME, the J-2S data are 
multiplied by a scale factor, ^ggj^g/j 2s' ^is sea-ling factor is also treated 

as a random variable. The distribution for the scaling factor is normal with 
mean equal to 1. 53 and a standard deviation of 0. 06. This information was 
supplied by S&E-ASTN-PPT, MSFC. 

In view of the above discussion the x and y components of a side load 
excitation acting on grid point 33 are expressed as 

F x = (ifih) (cose) (E SSME/J-2S> O- 1 ’ F i (t > (1 > 


and 

F y ^ = ( 1432 ) ( sin0 ) ( R SSME/J-2 s) F i ^ ^ 

where F. (t) is one of the 24 side load transients presented in Figures 7 through 

30; i.e., i= 1, 2, 3, . . . 24. The transient number, i , is also treated as 
a random variable uniformly distributed between 1 and 24. The factor 1. 1 
is included. to account for the decrease in peak values associated with digitiza- 
tion of the analog tape. 

It is observed that this description of the side load forcing function 
involves four variables F*, 9, -^gg^g/j 2s’ anc * F - (*")• Considering the 

ranges and various combinations of these variables, an infinitude of side load 
forcing functions can be defined. 


MONTE CARLO LOADS ANALYSIS 


There are at least two methods for determining dynamic loads and 
responses. One method is to calculate loads based upon discrete values of the 
input parameters selected by judgment from the entire range of possible input 
data. If judgment and experience are sufficient to identify the combination of 
parameters yielding maximum loads, the resulting maximum load may be 
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unduly conservative. For example, if the worst combination of the four 
variables involved in the SSME side loads analysis is selected as input, the 
resulting loads will exceed the design values. A second method is to utilize 
the Monte Carlo technique. This involves drawing samples of the input param- 
eters from their parent populations through sampling procedures governed by 
specified probability laws. 

For this analysis the input parameters which define the side load 
excitation are treated as random variables. These parameters are as follows: 

1. F*, the amplitude of the J-2S side load, lb rms per 1. 6 Hz band- 
width. 

2. 0, the direction of the side load excitation. 

3. R gsME/J 2S’ scaRn S f actor from J-2S data to SSME. 

4. F^ (t), the side load transient number, i, from the set of 24 
available. 

Values of the above variables were drawn to define 100 side load forcing 
functions. The procedure used herein for drawing the variables according to 
their probability laws is based upon a sequence of random numbers, p, , 

uniformly distributed in the unit interval 0 =£ x ^ 1. UNIVAC library subroutine 
RANDOM was used to generate 400 random numbers uniformly distributed 
between 0 and 1. Four of these are used to determine values of the input 
variables for each of the 100 cases. Having determined the random numbers, 
p. , to select samples of random variable, £ , with distribution F (x) one must 

solve the equation 


p = F (x) 


(3) 


for x as a function of p . The solution is 


x = 0 (p) . 


( 4 ) 
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Substitute the random numbers p. into the right side of equation (4) and 

obtain a set of numbers, t. , as 

1 

fi 


\ = <t> (P { ) 


( 5 ) 


The numbers are distributed according to the distribution F(x). 2 The input 

variables for this study defined using this random selection technique are 
presented in Table 5. The peak values of the resulting side load forcing 
functions are presented in Table 6. 

Using the 100 side load forcing functions generated by the random 
selections procedure, 100 transient response cases were run resulting in peak 
loads and accelerations for each case. The peak loads and acceleration data 
are combined statistically, resulting in a single set of loads to be used for 
strength assessment purposes. The combination of load is 


L = L ± KS, 


where 


L = statistical load, 


K *= number of standard deviations desired. 


S = sample standard deviation of the load 


N 

l (E - 




9 


2. Monte Carlo Principles and Neutron Transport Problems, by Jerome 
Spanier and Ely Gelbard. 
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L 


mean value of the N cases. 



y 


N = number of cases, (100 for this analysis) , 


and 


L. 

1 


peak value ( maximum or minimum) for case i. 


Since a finite sample size is used to form the statistical load (100 cases 
in this analysis) , confidence limits are applied to obtain an estimate of the upper 
limit on the true value of the statistical load. The confidence limits are 
determined by modifying the K value to assure that a load having K sigma prob- 
ability will not be exceeded. This new modified value of K is designated K . 3 

The load combination thus becomes 


L = L ± K S 
P 

where 



K = number of standard deviations (three for this analysis) , 


3. The expression for K is explained in Monte Carlo Dynamic Analysis for 

Lunar Module Loads, by Merchant and Sandy, Journal of Spacecraft and 
Rockets, Vol. 8, No. 1, January 1971. 
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Z = point on the normal density curve corresponding to the desired 
P confidence level, 

Zp = 1. 645 for 95 percent confidence level ( used in this analysis) 

and 

N = number of cases (100 in this analysis). 


Summarizing, the following steps were taken to determine the statistical 
loads and accelerations for the SSME side load excitation problem: 

1. 100 sets of randomly selected input parameters are chosen. 

2. 100 transient analyses are accomplished to obtain 100 cases of 

loads. 

3. Peak (maximum and minimum) loads are obtained and are assumed 
to have a Gaussian distribution. 

4 . A mean and standard deviation are computed. 

5. A statistical load is computed for strength analysis. This load is 
defined to be the upper confidence limit on the 3cr load corresponding to a 95 
percent confidence level. 


OUTPUT 


Since the three dynamic models of the engines are different, because of 
differing backup structure and orientations, they will necessarily yield different 
load responses. For expediency engine 3 was not considered; the justification 
for this is as follows: 

Previous analysis indicates that in general the loads tend to increase as 
the pendulum frequency increases. Engine 2 has the highest pendulum frequency 
and thus must be considered. Engine 1 being the upper engine is also considered 
for completeness and comparison. Thus one has loads for both a lower engine, 
engine 2, and the upper engine, engine 1. 
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The transient response analysis was accomplished using the 100 forcing 
functions generated as described above. The engines were in the start posi- 
tions, 4 which are identified in Figure 1. 

The response data for each engine was divided into two sets, output set 
A and output set B. The actuator loads and element loads are presented in out- 
put set A. Displacement and acceleration responses are in output set B. The 
coordinates ( or rows) of output sets A and B are identified in Tables 7 and 8, 
respectively. The statistical values of output set A and output set B for engine 
1 are presented in Tables 9 and 10, respectively. Similarly, output set A and 
output set B data for engine 2 are presented in Tables 11 and 12, respectively. 
The statistical values represent the upper confidence limit on the 3o value for 
a 95 percent confidence level as determined from the 100 case samples. The 
mean values and standard deviations are included for information only. For 
strength assessment the statistical values should be used. 


RESULTS 

Results of particular interest are extracted from the output tables and 
tabulated below. 

1. Actuator Load: 

Maximum = 69 980 lb (Compression), Table 9, row 2 

Minimum = 73 020 lb (Tension), Table 9, row 2 

2. Throat Bending Moment Load: 

Peak = 2 672 000 in. -lb, Table 9, row 43 

3. Gimbal Bearing Center Torsion Load: 

.Peak = 598 200 in. -lb, Table 11, row 13 


4. Rockwell Internal Letter, P&PS-GHC-73-029, April 26, 1973. 
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4. High Pressure Fuel Pump Acceleration: 


Translational: 

Peak = 25. 5 g, Table 10, row 22 
Rotational: 

Peak = 1. 23g/in., Table 10, row 26 

5. High Pressure Oxidizer Pump Acceleration 
Translational: 

Peak = 33. 06 g, Table 10, row 30 
Rotational: 

Peak = 1. 16 g/in., Table 10, row 34 

6. Nozzle Exit Displacement 

Peak = 2. 19 in.. Table 12, row 3 

The actuator incorporates a bypass which opens at an actuator load of 
90 000 lb. The 3a actuator load (73 020 lb) is well below the bypass value. It 
is interesting that for each engine only one case out of the 100 cases resulted in 
an actuator load exceeding the bypass value. 

The nozzle exit displacement of 2. 2 in. is well within the available rattle 
space of approximately 18 in. To demonstrate the transient character of the 
response data, time histories of actuator loads for a specific case are presented 
in Figures 31 and 32. These time histories are typical of those obtained for all 
cases examined. 
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La;;vliLa ACTUATOR 

Xq, Yq. Z q ORBITER COORDINATE AXES 

X, Y. Z ENGINE COORDINATE AXES 


ENGINE 

NULL POSITION 

START POSITION 

NUMBER 

PITCH 

YAW 

PITCH 

YAW 

i 

16 DEG 

ODEG 

18 DEG 

ODEG 

2&3 

10 DEG 

3.S DEG 

7 DEG 40 MIN 

S DEG 10 MIN 


NOTES: ENGINES CANTED OUTBOARD IN YAW 
ENGINES CANTED BELL UP IN PITCH 

Figure 1. SSME-orbiter configuration. 
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K BUS • ORBITER BACKUP STRUCTURE STIFFNESS MATRIX 


K d - DUCT STIFFNESS MATRICES 


Figure 3. SSME structural model. 
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SLOPE = 0.01 cis/in.-lb 
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Figure 5. Servoactuator schematic. 
















Figure 6. Distribution of J-2S rms amplitude for 1.6 Hz bandwidth. 
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Side Load Forcing Functions are generated using from this set 
of 24 transients (Figs. 7 through 30) as follows: 

F.* 

F j ^ = 1432 ( R SSME/J2 -s) F i ^ ;i= 1,2,3 24 

j = 3» » • # ® 
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Figure 7. Side load transient no. 1. 
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Figure 9. Side load transient no. 3. 
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Figure 10. Side load transient no. 4. 


21 


ii 


ii ..M 


II 11 .11 M 11 Ii .11 11 


u.. u 


ii ii U 11 L L . L 


AMPLITUDE (lb) 



Figure 11. Side load transient no. 5. 
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Figure 17. Side load transient no. 11. 
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Figure 18. Side load transient no. 12. 
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Figure 19. Side load transient no. 13. 
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Figure 21. Side load transient no. 15. 
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Figure 22. Side load transient no. 16. 
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Figure 23. Side load transient no. 17. 
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Figure 24. Side load transient no. 18. 
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Figure 25. Side load transient no. 19. 
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Figure 26. Side load transient no. 20. 
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Figure 27. Side load transient no. 21. 
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Figure 29. Side load transient no. 23. 
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Figure 30 . Side load transient no. 24 . 
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TIME (SEC) 

Figure 31. Actuator no. 1 load time history engine 2, start position, 
excitation case 1 peak side load excitation = 13 613.3 lb. 
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Figure 32. Actuator no. 2 load time history engine 2, start position, 
excitation case 1 peak side load excitation = 13 613.3 lb. 
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TABLE 1. DEGREE OF FREEDOM TABLE 


— 

NASTRAN 

Node 

Number 


Degree of Freedom 

Node Location (in.) 

NASTRAN 

Coord. 

System 

Node 

Description 

X 

Y 

Z 

RX 

RY 

RZ 

X 

Y 

Z 

1 

1 

2 

3 

4 

5 

6 

0.0 

0.0 

0.0 

4 

Gimbal Point 

2 

7 

8 

9 

10 

11 

12 

- 15. 5 

0.0 

0.0 

3 

Main Beam/ Pumps 

3 

13 

14 

15 

16 

17 

18 

-27.9 

0.0 

0.0 

3 

Thrust Chamber C. L. 

4 

19 

20 

21 

22 

23 

24 

-39.9 

0.0 

0.0 

3 

Throat C. L. 

5 

25 

26 

27 

28 

29 

30 

-29.4 

-19.5 

-19.5 

3 

HPFTP 

6 

31 

32 

33 

34 

35 

36 

-28.6 

18.4 

18.4 

3 

HPOTP 

7 

37 

38 

39 




10.0 

-31.5 

0.0 

4 

Actuator l/Orbiter 

8 

40 

41 

42 




10.0 

0 . 0 

-31.5 

4 

Actuator 2/ Orbiter 

9 

43 

44 

45 

46 

47 

48 

-36.0 

-25. 0 

0 . 0 

3 

Actuator l/Struts 

10 

49 

50 

51 

52 

53 

54 

-36.0 

0.0 

-25. 0 

3 

Actuator 2/Struts 

11 

55 

56 

57 

58 

59 

60 

-45.75 

0.0 

0.0 

3 

Main Beam/ Struts 

12 

61 

62 

63 

64 

65 

66 

-24. 88 

-10. 72 

7.39 

3 

Strut 

13 

67 

68 

69 

70 

71 

72 

-24. 88 

-10. 72 

-7.39 

3 

Strut 

14 

73 

74 

75 

76 

77 

78 

-24. 88 

7.39 

-10. 72 

3 

Strut 

15 

79 

80 

81 

82 

83 

84 

-24.88 

7.39 

-10.72 

3 

Strut 

16 

85 

86 

87 

88 

89 

90 

-45.75 

10.20 

-5.89 

3 

Strut 

17 

91 

92 

93 

94 

95 

96 

-41. 52 

-12. 12 

7. 04 

3 

Strut 

18 

97 

98 

99 

100 

101 

102 

-41.52 

-12. 12 

-7. 04 

3 

Strut 

19 

103 

104 

105 

106 

107 

108 

-45. 75 

-5.89 

-10.20 

3 

Strut 

20 

109 

110 

111 

112 

113 

114 

-45. 75 

-10.20 

-5.89 

3 

Strut 

21 

115 

116 

117 

118 

119 

120 

-41. 52 

7. 04 

-12. 12 

3 

Strut 

22 

121 

122 

123 

124 

125 

126 

-41. 52 

-7. 04 

-12. 12 

3 

Strut 

23 

127 

128 

129 

130 

131 

132 

-45. 75 

-5.89 

10.2 

3 

Strut 

24 

133 

134 

135 

136 

137 

138 

-54. 0 

0 . 0 

0.0 

3 

Bell C. L. 

25 

139 

140 

141 

142 

143 

144 

-63. 13 

0.0 

0.0 

3 

Bell C. L. 

26 

145 

146 

147 

148 

149 

150 

-74. 19 

0.0 

0 . 0 

3 

Bell C. L. 

27 

151 

152 

153 

154 

155 

156 

-84.99 

0 . 0 

0.0 

3 

Bell C. L. 

28 

157 

158 

159 

160 

161 

162 

-95. 79 

0.0 

0.0 

3 

Bell C. L. 

29 

163 

164 

165 

166 

167 

168 

-106.59 

0 . 0 

0.0 

3 1 

Bell C. L. 

30 

169 

170 

171 

172 

173 

174 

-117.59 

0.0 

0.0 

3 

Bell C. L. 

31 

175 

176 

177 

178 

179 

180 

-129.59 

0.0 

0 . 0 

3 

Bell C. L. 

32 

181 

182 

183 

184 

185 

186 

-142.49 

0.0 

0.0 

3 

Bell C. L. 

33 

187 

188 

189 

190 

191 

192 

-155.59 

0.0 

0.0 

3 

Bell C. L. at Exit 

34 

193 

194 

195 

196 

197 

198 

-20.0 

-4.55 

9. 11 

3 

Strut 

35 

199 

200 

201 

202 

203 

204 

-20.0 

-4.55 

-9. 11 

3 

Strut 

36 

205 

206 

207 

208 

209 

210 

-20.0 

9. 11 

-4.55 

3 

Strut 

37 

211 

212 

213 

214 

215 

216 

-20.0 

-9. 11 

-4.55 

3 

Strut 

38 

217 

218 

219 

220 

221 

222 

-34.8 

-12.0 

26.0 

3 

Controller C.G. 

39 

223 

224 

225 

226 

227 

228 

-20.0 

0.0 

0.0 

3 

Main Beam/Struts 

40 

229 

230 

231 

232 

233 

234 

1.5 

21. 52 

-21. 26 

4 

LPOTP 

41 

235 

236 

237 

238 

239 

240 

-1.71 

-21. 0 

21.0 

4 

LPFTP 

43 

241 

242 

243 

244 

245 

246 

10.0 

11.0 

12.54 

1 3 

LPFTP Duct 

44 

247 

248 

249 

250 

251 

252 

-7.56 

-16.3 

-16.3 

3 

LPFTP Duct 

45 

253 

254 

255 

256 

257 

258 

-8.74 

16.68 

16.09 

3 

LPFTP Duct 


NOTES: Use Right Hand Rule 

Grid Points are identified in Figure 2. 
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TABLE 2. IDENTIFICATION OF SSME DUCTS 


•NASTRAN 

Duct 

No. 

Duct 

Discription 

Nodes 
Connected 
By Duct 

1 

LPFTP Turbine Drive 

39, 41 

2 

LPFTP Discharge 

41, 5 

3 

LPOTP Turbine Drive 

40, 6 

4 

LPOTP Discharge 

40, 6 

5 

LPFTP Turbine Discharge 

11, 41, 43, 44, 45 


TABLE 3. IDENTIFICATION OF STRUCTURAL ELEMENTS 
SELECTED FOR LOADS OUTPUT 


Element 

Number 

Node 
At End A 

Node 
At End B 

Nodes Defining 
Plane 1 

Type 

Element 

1 

1 

2 

1 

2 

7 

Bar 

2 

2 

39 

2 

39 

7 

Bar 

3 

3 

4 

3 

4 

7 

Bar 

4 

4 

11 

4 

11 

7 

Bar 

9 

5 

2 

5 

2 

3 

Bar 

10 

6 

2 

6 

2 

3 

Bar 

11 

39 

3 

39 

3 

7 

Bar 

12 

11 

24 

11 

24 

7 

Bar 

13 

24 

25 

24 

25 

7 

Bar 


NOTE: Element coordinate system defined in Figure 4. 
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TABLE 4. SERVOACTUATOR PARAMETERS 


Symbol 

Units 

— • 

Value 

Description 

K 

V 

cis(in.-lb) 

1228 


Servovalve flow gain 

G 

V 

W 

V 

rad/sec 

w 

220 

-1 

\ "1 
( — L] s+ 1 
V w / 

\ v/ J 

Servovalve dynamics 
Servovalve natural frequency 


- 

0.7 


Servovalve damping coefficient 

A 

in. 2 

24.83 


Actuator piston area 

K 

o 

Ib/in. 

902 800 


Oil stiffness 

M 

P 

( lb-sec 2 )in. 

0. 11 


Actuator piston equivalent mass 

K 

P 

cis/psi 

0. 15 


Pressure feedback gain 



T S 
P 


Pressure feedback shaping network 
Pressure feedback time constant 
Position feedback gain 

G 

P 

T 

P 

H 

sec 

( in.-lb)/in. 

T S + 1 
P 

0. 125 
0. 364 


f 

lb 

- 


External disturbance load 

f a 

lb 

- 


Actuator developed force 

f b 

lb 

90 000 


Bypass actuating force 

X. 

1 

in. 

— Q 

A ^0 ^Adt 


Ideal piston position 

X 

p 

in. 

- 


Actuator piston position 

Q v 

cis 

- 


Servovalve flow 

Q dpf 

cis 

- 


Dynamic pressure feedback flow 

q b 

cis 

- 


Load bypass flow 


cis 

- 


Flow to actuator piston 

E 

in.- lb 

- 


Actuating error signal 

T 

c 

in. -lb 

- 


Command torque 
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TABLE 5. RANDOMLY SELECTED INPUT DATA 



Transient 

J-2S Side Load Amplitude 

SSME Scaling 

Direction of 

Case 

Number 

[lb rms/SQKT (1.6 Hz)] 

Factor 

Transient (deg) 

1 

22 

940. 90 

u 

19:03 

2 

3 

337.13 


72.89 

3 

19 

186.18 


1/73 

4 

6 

461.81 


84.95 

5 

13 

463.43 


277.08 






6 

7 

St » SI 

1. 537 

336. 69 

7 

16 

104. 58 

1.468 

302.16 

8 

16 

778. 26 

1. 571 

293.18 

9 

1 

740. 72 

1.453 

282.41 

10 

16 

151.04 

1.500 

345.91 

11 

12 

106.60 

1.556 

326. 82 

12 

1 

600. 96 

1.483 

231.09 

13 

15 

292.06 

1.575 

19.38 

14 

16 

872.84 

1.589 

204.27 

15 

16 

520. 91 

1.485 

206.75 

16 

22 

134.04 

1.477 

291.12 

17 

9 

554.28 

1.454 

162.26 

18 

21 

304.70 

1.553 

134.52 

19 

3 

976.57 

1.399 

206.06 

20 

3 

594.23 

1.426 

157.88 

21 

10 

165.28 

1.496 

61.55 

22 

16 

438.45 

1.461 

84.68 

23 

16 

107.19 

1.419 

148.15 

24 

8 

402. 82 

1. 536 

214.33 

25 

24 

418. 23 

1.469 

50.16 

26 

23 

154.87 

1.451 

54.24 

27 

21 

279. 96 

1.487 

282.02 

28 

15 

460. 55 

1.531 

68.09 

29 

8 

375.64 

1.505 

149. 71 

30 

9 

345.75 

1.536 

60.48 

31 

4 

664.51 

1.562 

348.47 

32 

22 

450.54 

1.531 

317.65 

33 

18 

114.74 

1.513 

216.75 

34 

1 

458.54 

1.489 

104.65 

35 

19 

415.76 

1.575 

86.37 

36 

15 

305. 61 

1. 523 

35.86 

37 

18 

541.42 

1.652 

89.02 

38 

24 

263. 59 

1.521 

138.29 

39 

24 

376. 66 

1.484 

98. 91 

40 

3 

370.22 

1.497 

149.95 

41 

23 

703. 09 

1.415 

91.72 

42 

19 

452.00 

1.523 

29.40 

43 

i9 

100.98 

1.515 

199.83 

44 

7 

308.64 

1.556 

180.31 

45 

4 

589.17 

1.597 

103.72 

46 

23 

371.57 

1.514 

248.78 

47 

8 

449.48 

1.518 

279.70 

48 

22 

599. 09 

1.539 

354.18 

49 

11 

298. 55 

1. 515 

63. 83 

50 

15 

895.35 

1.577 

36.19 
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TABLE 5. (Concluded) 


Case 

Transient 

Number 

J-2S Side Load Amplitude 
[ lb rms/SQRT (1.6 Hz)] 

SSME Scaling 
Factor 

Direction of 
Transient (deg) 

51 

16 

203.35 

1.414 

85. 82 

52 

19 

84.97 

1.516 

40.76 

53 

15 

518.47 

1.614 

114. 89 

54 

2 

728. 92 

1.489 

311.79 

55 

19 

144. 89 

1.549 

317. 77 

56 

3 

370. 75 

1.553 

209.77 

57 

23 

532.00 

1.526 

150.22 

58 

1 

169.76 

1.589 

70.15 

59 

9 

152.56 

1.563 

252.67 

60 

11 

604. 77 

1. 535 

137.32 

61 

17 

113.44 

1.522 

332.97 

62 

18 

433. 57 

1. 504 

354.49 

63 

10 

310.62 

1.552 

43.08 

64 

13 

126.22 

1.570 

116.17 

65 

20 

463.92 

1.672 

237. 58 

66 

10 

406.72 

1.489 

195.07 

67 

11 

1578.34 

1.528 

134.12 

68 

1 

509.06 

1.582 

85.05 

69 

17 

282.84 

1.657 

134.79 

70 

18 

413.87 

1.605 

189.57 

71 

14 

336. 89 

1.581 

112.48 

72 

13 

306. 77 

1.448 

351.30 

73 

2 

293.28 

1.498 

290.42 

74 

17 

221. 88 

1.484 

45.90 

75 

14 

466.38 

1. 568 

126.19 

76 

14 

426.49 

1.490 

324.48 

77 

15 

125.42 

1.525 

265.62 

78 

15 

156.58 

1. 506 

306. 82 

79 

12 

302.42 

1.519 

198.07 

80 

24 

437.73 

1.580 

67.60 

81 

17 

202. 84 

1. 536 

314.14 

82 

15 

477.19 

1.680 

85. 61 

83 

20 

201.82 

1. 554 

0.53 

84 

24 

198. 89 

1.662 

337. 56 

85 

4 

610.48 

1.578 

176.72 

86 

1 

397.79 

1.589 

80.03 

87 

20 

374.13 

1.614 

264.67 

88 

12 

384.98 

1.494 

251.20 

89 

5 

88. 59 

1.518 

337.99 

90 

9 

.339. 63 

1.483 

164.78 

91 

21 

545.05 

1.475 

201.36 

92 

17 

528.05 

1.570 

218. 79 

93 

16 . 

390. 83 

1.471 

56.33 

94 

10 

568. 83 

1.517 

126.21 

95 

24 

134.86 

1.549 

117.29 

96 

22 

483.47 

1.577 

112.23 

97 

20 

511.01 

1.441 

201.66 

98 

19 

385.77 

1.504 

88. 73 

99 

4 

488.68 

1.393 

93.59 

100 

13 

522.69 

1.596 

52.67 
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TABLE 6. PEAK FORCE OF SIDE LOAD TRANSIENTS 


Case 

Transient 

Number 

Peak Force 
(ib) 

Case 

Transient 

Number 

Peak Force 
(ib) 

1 

22 

13613.3 

51 

16 

24 3 6.9 

2 

3 

4278.9 

52 

19 

1216.3 

3 

19 

2567.4 

53 

15 

9750.6 

4 

6 

7028. 9 

54 

2 

10601.6 

5 

13 

7142.1 

55 

19 

2118.8 

6 

7 

1450.7 

56 

3 

4957.2 

7 

16 

1301.3 

57 

23 

8484.3 

8 

16 

10359.3 

58 

1 

2920. 5 

9 

1 

11650.1 

59 

9 

2801.6 

10 

16 

1920.3 

60 

11 

10215.3 

11 

12 

1719.2 

61 

17 

2105.7 

12 

1 

9644.3 

62 

18 

6426.8 

13 

15 

5360.0 

63 

10 

5102.1 

14 

16 

11750.0 

64 

13 

2026.3 

15 

16 

6555.0 

65 

20 

8105.8 

16 

22 

1961.0 

66 

10 

6410.3 

17 

9 

10348.7 

67 

11 

26546.8 

18 

21 

4531.5 

68 

i 

8715.9 

19 

3 

11771.2 

69 

17 

5716.3 

20 

3 

7300.4 

70 

18 

6544.2 

21 

10 

2617.2 

71 

14 

4784.1 

22 

16 

5429.1 

72 

13 

4541.5 

23 

16 

1289.4 

73 

2 

4287.7 

24 

8 

5555.5 

74 

17 

4015.8 

25 

24 

7641.3 

75 

14 

6568.9 

26 

23 

2349.2 

76 

14 

5708.4 

27 

21 

3985.2 

77 

15 

2228. 6 

28 

15 

8218.1 

78 

15 

2747.6 

29 

8 

5075.8 

79 

12 

4760.2 

30 

9 

6820.0 

80 

24 

8598. 9 

31 

4 

11374.7 

81 

17 

3799.1 

32 

22 

6830.3 

82 

15 

9341.3 

33 

18 

1710.9 

83 

20 

3278.0 

34 

1 

7387.9 

84 

24 

4110.6 

35 

19 

6180.7 

85 

4 

10560.8 

36 

15 

5421.6 

86 

1 

6840.8 

37 

18 

8814.9 

87 

20 

6313.5 

38 

24 

4983.0 

88 

12 

5961.9 

39 

24 

6947. 8 

89 

5 

1579.3 

40 

3 

4773.7 

90 

9 

64 65.8 

41 

23 

10397.9 

91 

21 

7695.2 

42 

19 

64 98.0 

92 

17 

10111.6 

43 

19 

1444.5 

93 

16 

4871.9 

44 

7 

4777.3 

94 

10 

9135.9 

45 

4 

10310.5 

95 

24 

2596.6 

46 

23 

5879.5 

96 

22 

7547.6 

47 

8 

6126.6 

97 

20 

7697.0 

48 

22 

9128.5 

98 

19 

5477.1 

49 

11 

4978.0 

99 

4 

7461.5 

50 

15 

16444.3 

100 

13 

8526.4 
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TABLE 7. IDENTIFICATION OF OUTPUT SET A 


Row 

Name 

Element 

End 

Plane 

Node 

Description 

1 

Act. 1 





Actuator 1 Load 

2 

Act. 2 





Actuator 2 Load 

3 

Ml A1 

1 

A 

1 

1 

Gimbal Point, Plane 1 Moment 

4 

Ml A2 

1 

A 

2 

1 

Gimbal Point, Plane 2 

5 

Ml A3 

1 

A 

3 

1 

Gimbal Point, Resultant Moment 

6 

Ml B1 

1 

B 

1 

2 



7 

Ml B2 

1 

B 

2 

2 



8 

Ml B3 

1 

B 

3 

2 



9 

SI 1 

1 


1 


Gimbal Point, Plane 1 Shear 

10 

SI 2 

1 


2 


Gimbal Point, Plane 2 Shear 

11 

SI 3 

1 


3 


Gimbal Point, Resultant Shear 

12 

Al 

1 




Gimbal Point, Axial Load 

13 

T1 

1 




Gimbal Point, Torsion (-X Side) 

14 

M2 Al 

2 

A 

1 

2 

Engine C. L 

. Loads 

15 

M2 A2 

2 

A 

2 

2 



16 

M2 A3 

2 

A 

3 

2 



17 

M2 B1 

2 

B 

1 

39 



18 

M2 B2 

2 

B 

2 

39 



19 

M2 B3 

2 

B 

3 

39 



20 

S2 1 

2 


1 




21 

S2 2 

2 


2 




22 

S2 3 

2 


3 




23 

A2 

2 






24 

T2 

2 






25 

Mil Al 

11 

A 

1 

39 



26 

Mil A2 

11 

A 

2 

39 



27 

Mil A3 

11 

A 

3 

39 



28 

Sll 1 

11 


1 




29 

Sll 2 

11 


2 




30 

Sll 3 

11 


3 




31 

All 

11 






32 

Til 

11 






33 

M3 Al 

3 

A 

1 

3 



34 

M3 A2 

3 

A 

2 

3 



35 

M3 A3 

3 

A 

3 

3 



36 

S3 1 

3 


1 


I 


37 

S3 2 

3 


2 


\ 


38 

S3 3 

3 


3 


Engine C. L 

. Loads 
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TABLE 7. (Continued) 


Row 

Name 

Element 

End 

Plane 

Node 

Description 

39 

A3 

3 




Engine C. L. Loads 

40 

T3 

3 




Engine C. L. Loads 

41 

M4 A1 

4 

A 

1 

4 

Throat, 

Plane 1 Moment 

42 

M4 A2 

4 

A 

2 

4 

Throat, 

Plane 2 Moment 

43 

M4 A3 

4 

A 

3 

4 

Throat, 

Resultant Moment 

44 

M4 B1 

4 

B 

1 

11 

Engine C. L. Loads 

45 

M4 B2 

4 

B 

2 

11 



46 

M4 B3 

4 

B 

3 

11 



47 

S4 1 

4 


1 




48 

S4 2 

4 


2 




49 

S4 3 

4 


3 




50 

A4 

4 






51 

T4 

4 






52 

M12A1 

12 

A 

1 

11 



53 

M12A2 

12 

A 

2 

11 



54 

M12A3 

12 

A 

3 

11 



55 

S12 1 

12 


1 




56 

S12 2 

12 


2 




57 

S12 3 

12 


3 




58 

A 12 

12 






59 

T12 

12 




Engine C. L. Loads 

60 

M13A1 

13 

A 

1 

24 

Bell C. L. Loads 

61 

M13A2 

13 

A 

2 

24 



62 

M13A3 

13 

A 

3 

24 



63 

M13B1 

13 

B 

1 

25 



64 

M13B2 

13 

B 

2 

25 



65 

M13B3 

13 

B ! 

3 

25 



66 

S13 1 

13 


1 




67 

S13 2 

13 


2 




68 

S13 3 

13 


3 




69 

A 13 

13 




1 

f 

70 

T13 

13 




Bell C. L. Loads 

71 

M9 B1 

9 

B 

1 

2 

HPFTP Beam Loads 

72 

M9 B2 

9 

B 

2 

2 



73 

M9 B3 

9 

B 

3 

2 



74 

S9 1 

9 


1 


1 

1 

75 

S9 2 

9 


2 


1 

r 

76 

S9 3 

9 


3 


HPFTP Beam Loads 
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TABLE 7. (Concluded) 


Row 

Name 

Element 

End 

Plane 

Node 

Description 

77 

A9 

9 




HPFTP Beam Loads 

78 

T9 

9 




HPFTP Beam Loads 

79 

M10B1 

10 

B 

1 

2 

HPOTP Beam Loads 

80 

M10B2 

10 

B 

2 

2 



81 

M10 B3 

10 

B 

3 

2 



82 

S10 1 

10 


1 




83 

S10 2 

10 


2 




84 

S10 3 

10 


3 




85 

A10 

10 






86 

T10 

10 




HPOTP Beam Loads 


NOTES: Prefixes: M = Bending Moment ( in. -lb) 

S = Shear Load ( lb) 

A = Axial Load ( lb) 

T = Torsional Load (in.-lb) 

Plane 3 loads are resultant of loads for plane 1 and plane 2; for 
example: M10B3 = SQRT( M10B1**2 + M10B2**2) 
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TABLE 8. 


Row 

Name 

1 

DY33 

2 

DZ33 

3 

DYZ33 

4 

AX1 

5 

AY1 

6 

AZ1 

7 

AYZ1 

8 

RAX1 

9 

RAY1 

10 

RAZ1 

11 

AX2 

12 

AY2 

13 

AZ2 

14 

AYZ2 

15 

RAX2 

16 

RAY2 

17 

RAZ2 

18 

RAYZ2 

19 

AX5 

20 

AY5 

21 

AZ5 

22 

AXYZ5 

23 

RAX5 

24 

RAY5 

25 

RAZ5 

26 

RAXYZ5 

27 

AX6 

28 

AY6 

29 

AZ6 

30 

AXYZ6 

31 

RAX6 

32 

RAY6 

33 

RAZ6 

34 

RAXYZ6 
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Node 


Description 


33 

33 

33 

1 

1 

1 

1 

1 

1 

1 

2 

2 

2 

2 

2 

2 

2 

2 

5 

5 

5 

5 

5 

5 

5 

5 

6 
6 
6 
6 
6 
6 
6 
6 


Bell Displacement at Exit 

♦ 

Bell Displacement at Exit 
Gimbal Point Accelerations 


Gimbal Point Accelerations 


HPFTP Accelerations 


» 

HPFTP Accelerations 
HPOTP Accelerations 


t 

HPOTP Accelerations 









TABLE 8. (Concluded) 


Row 

Name 

Node 

Description 

35 

AX 11 

11 


36 

AY11 

11 


37 

AZ 11 

11 


38 

AYZ11 

11 


39 

RAX11 

11 


40 

RAY11 

11 


41 

RAZ11 

11 


42 

RAYZ 11 



43 

AX38 

38 

Controller Accelerations 

44 

AY38 

38 



45 

AZ38 

38 



46 

AXYZ38 

38 



47 

RAX38 

38 

>• a 


48 

RAY38 

38 

T 

49 

RAZ38 

38 

Controller Accelerations 


NOTES: Prefixes: DY 
DZ 
DYZ 
AX 
AY 
AZ 
AYZ 
AXYZ 
RAX 
RAY 
RAZ 
RAYZ 
RAXYZ 


Displacement in Y-Direction 

Displacement in Z-Direction 

Root Sum Square of Displacements, DY, DZ 

Acceleration in X-Direction 

Acceleration in Y-Direction 

Acceleration in Z-Direction 

Root Sum Square of Accelerations, AY, AZ 

Root Sum Square of Accelerations, AX, AY, AZ 

Rotational Acceleration About X-Axis 

Rotational Acceleration About Y-Axis 

Rotational Acceleration About Z-Axis 

Root Sum Square of Rotational Acceleration, RAY, RAZ. 

Root Sum Square of Rotational Acceleration, RAX, RAY, RAZ. 


Displacements in inches, Translational Accelerations in g, Rotational Accelerations in g/inch 
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TABLE 9. STATISTICAL VALUES — OUTPUT SET A, ENGINE 1 — START POSITION 
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TABLE 9. (Concluded) 
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TABLE 10. STATISTICAL VALUES - OUTPUT SET B, ENGINE 1 - START POSITION 
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I ZOO OOOOOOOOOOOOOOOOCOOC'COCOCCOC'COOOCCOOOOCOCO cccc 

IT. I e O ♦ 4 4 4 4 1 14 4 4 4 4 I I 44 I 4 4 4 I I 4 4 I 4 4 4 I I 4 4 I 4 4 4 I I I I 14 4 4 I I I I 

z lorHCDificw^coirtfHCvco'w^ccoc'W^ccHXtrwcvirtrc^cccctv'Hcif-cwwir^h^O' •- cc ir in 

r i«kMfccir^^if»rific(\co»f'f'COffcoi r i< , i f tcf< , if.ircirh\£h<\^^r^^»rK.fth * — 

s i c«z^ccccK:ir^<\jir*Mcc^crtrircv.caccaK:r^^cc’<vvc — vcv£hv£cc^ccjir»r>ccc * * r~ cr^cc 
»- i z^K 2 K;o^ 4 rc^^HC\iOJ^-iwoj<v*-!C\j^— NC-^^tvm<V(v^H^\£cv<\irc<\.«-i^ 4 *-^ — ^cc^accrorc— («~i — a a: 

Z I « > 

►* I K U 
5 I IT C 
I 

S 1 

I 

| o O HHrHlOHOOWHHlOHOO^H 4 H(\JHOOiH 4 H(\lHCOiHH 4 lOHOHUnHH-( (\IHSH 

I C C CCCeoeeceeCCCCCCCCCCCCCCCecCCCCCCCCCCCCCCC cccc 

I , + 4 444 I 144444 I | 4 4 I 4 4 4 I I 4 4 I 4 4 4 I |44I444II4||444 114 1 

i ! »rtrcc\imcrsCcvmw%cir\C^O'cr-irc^f^»r;»r:^HirK:csCiraccrr^'*-*-it^r-^ 0 ‘r^ir. irai^ r- «- a •- 
i L.irwc^c^Nvor**iDwccc^ccccirirccK:r-fr^p-cac(vo'K:i^(Oir^<cvcccc<\icccrK:o s ^ X tv tr 
i 2 rc<cc\Cir»"*-c^ftMr^fficir. atfiftfdcsircvtfiro'ffcff^r^hMra %r — c a ir - x cvr«-c 

| < J^^'CCM^inK)fOK)^eV!C\JK)CJ»O^H«H^U , ;CVrO^r 0 <V(VJK;(C«VlC«-ifOC\lOJd-<V;«- 4 CV*-CL-H 0 '^£)d'<VCV HWHff 

! s > r r * r r r * r r r i # r r * i* r i* * r r r * r r r * / r r * r r r * r r r * r r »* * r r r * i" r r 

i 

i 


I _J HHi 4 HMN(\IC^C^(Ma'C\,CCCC(V«f\J\ l CCC^f\){V(V>\CCCwH^-JH».C C C <\J — C\ ^“ c r 

I < CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCXCC CCwC 

I «J 44444444444444444444444444444444444444 1444444 41+4 

I ►- r^c^vc<ccctccvwvc^r^irNvci^^O‘ 0 > h‘ — tfavCMrff^ffvcirirmnrca<rchWhtfhh c — vC»r 

I u HvC\CO'(VO'J‘C»-c*i^H^?’ l r! 0 ‘O i .:'rir ^ r ■o c m c ^ * jd ir c x v: r c* cc r- ro r"- r- r- ,r cv x r- 

i cr tr. ^^-r-ccctr. ccv^cr^-ior-^^f^r^^tc^ir^acccv^-c^-rr k. rcrc-i^cvrrccr-rccrvccc^ c«-r*-K: 

i^^.^^cvcc^^cv — — : r, •cc\;<r^cv^x^^rG— — r- x x <g 

I H 3 

I « _I 
I K «S 


I 

I 

I 

IT icC 4444444 < f 474444744444447 X 444444 f 4444444 r 4 i ’4444 + 7+7 

z itt + ira^hiir^iTMrK; ^irx cvrrc^irh-tvtvjit.xffx-o'cvirsc^^aaxa-cK'.irxxh'tv c it k. it 

r i cKirrc^op^cvrcccxrrcr'crctxvrcvcr-a a , r~^a'Xv£d > r-r-r-crp |r :irir^d-^^<xa'f^vra c r- — ^ c 

z i o<f\j'j'c\jr* 3 'rom^\ 0 c\ix\ox x> ^ ou x inn X) -o 2 MA 1 -ts o r» xx 3 m 

t- 1 z+-rc<vj-^»ir^ir. m ^ ? h a mt, cv — + ^ if mt cv* »- a- ^ vf ^ a x ^ ^ o: c: ^ h h n h r ^ - -o ~ r- 

> 1 < > * 

<3 I K It 

St ts. c 

I 

I 

I 

I 

I ^ _ ^ 

I cococccocceccc cccc cccc ccccccccccocccccoc cccc C cccc 

| 44444441444444 14444444 I 4444+44 I 4444444 I 4 14444 4141 

i r- cv — \C cc x cc *r <v tr ^ ~ cv r- -- c d c r- c c - w <\ ■ & <r x rv s ? r- j* vr ^ j. x >c c c cv c c y r. «- r r 

I iLO^CtlO^^irc^MJ'^vClOWiTWHr.NO^'Xir^H.O^iO^-HNXWA'N^WNd-CircNiO ^ x u- c 

i ZDotnH^xtMiH^wxo'sin^^^d'vflMir. ^c\!K)H^ct^j , Hio^NxrtciriHN 1 nirovc x-n-^r- 
i < _t ^ ^ r- tv 3 t id x H-; cv cv k. ^ re ^ ^ ^ cv ro x x: cv cv rc- x rr rr O fC (\. ^’ ro cv ^ rv to cu x rv rw a~ cm x 

i u 

I s > 

I 

I 


i r vC at 

to tv it ki vC cm r x cr c 

Lui n m tn ^ ^ ^ ^ wcjcutvx ,nj in m lc >- n ^ x > h ^ n h ^ h h \ i x t- \ k. k- ; ^ 

Z I K]mNHH^NX 4 NW(\)(\JNX>-Ny^^ir;>-XVNX>r>CvC>-X)-NXHHHNX>.N>-mOr; >- X > f\. 

<r i >-rsi>-x>“isfV«r<<x>-NJ>-<i<t* , s<tx>-fs;x<i<i<<xvrsj>c<a:<i<i<x>-rs.>^<<<r<tx>-fV x «■ < < 


3 : i 4 «r. 4 tf:^N«ffc + ft.r. difXMtffc^K(ri(r^Mtffc^aK! 4 tf:«fhtfC , c^<viri«“ vtf^xc 

Cl ^^^^^^^^-4^rv.'f\JCVCVCVC\ ( fVCVCVOjrOir)PO’OfOrOfOrO?0»Oj--4^ ^ cfcrj- 3- 

cr i 
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TABLE 11. STATISTICAL VALUES - OUTPUT SET A, ENGINE 2 - START POSITION 


i . ' l i 

I « cccccccccccc©ccccccccccc©ccccccc,ccccccccccccc 

I O ++ + + I + + + + + + + + + + + + + + + + + + + + + ♦+ ++ + ♦’-♦♦ + +.+ + + + + + + +’+ 

i »- r^iOKiirvcar^ro^f^ttccojfncvnir^HCcv^ztacvcNjcirirr^xCcrccir^^^iriccecNDcvccitw 

i k aor-cc^v£)-HLrxCCT'^^ffir)crr^r^^^cvcccf^<v<vcvc^r-0'fo^(Tir<rcv*^'0»o^irr-ir«\j^ 

i t r. ^N^xccjcamoj^iriro'd’irtrd'vcccar^tvjcy'^'ra'irccwf^rrircco'^ca^cirirr-O'irf^ 

I fr- “ 

I < -II I I I I I I I I I I I I I I I I I I 1,1 I I I I I I I I I I I I I I I I I I I I I I I I 

I K < 
it n> 

l ; ; 

i : 

I 


ire*-3'sC'C<NJirin.^ir. sC\C\ct^sC'£t r )zt a- ~&\C\c*ct r izta~* ztyCsCsCi f )**-4&\C\CyCtn\C\c 
(©COOOCOCOOCOO.OOOOOCOOOCCOOCOCCCOCOOCCOCCCOCO 

+ + I I + + +•+ + + + + +'+ + ♦ + + + + + + + + + + ++ + + + .+ + + + + + + + + + + + + 

x tv ctrcu:^rf*-*“<v^ar^sro s u’c\ c^rr^^irtvct^cciriricvcrr^cff r r. ir — <\ r* cr 

<V,vC«AQ'Kir , «(J'K', O'ff^ICK', fC(V'0 O' a C C If, ^ (V a j ec c-wir^ — i^r-ctr 

(\,^(v^pr^HfOK:cv^(v^Hro^»*^rocN.irK'> 


~«U , >sOX)(\JLnm©ir)*CvO\DCVJX)X><\ja':*O3-X>X>vO<\J;±3'O^\0'£\C<\Jd‘.d’O^sGX)»U<\J\OvO 
‘cccceccecoccceecccecccecccccc'ccccccceccccc 
I I + ♦ + + + + + + + ♦+ + ♦+ ♦ + + + ♦ + ♦ + ++ + •♦.♦.« + + + + + + + + ♦+ + + + 
trr^acf^trccf^aa»rr*cK:r^c\.af^r^cccc<v^^H^ct\£C4rh'»r©c:<r\Ccv : (vr-^zt 
'r^in<\jintrccK)s£)^cvcc»occmf^c^-^ccKCh-h*K)oco'cco'r^^vCiro'0'i^cCif^O'a'c<cc 
:^^-cv^^trdifi^cc^ir^irtf^flrcv\Cffcc^r^^\ra^«actd-iriri^crd'i^cff\cor 
ih'\£ i c\j<v»H^H^-'r^*-t«rHcr^co^tc}'^irLr»^r-roir. Lr^-fOro^ir^j-ir rororo-Hcr^-t^d-rr^-^- 


t 


i _J LfircccvCvCMT, if^ifvTNNNNNMT, ir in^f^i^p^f^irirx ir. vcr^r^r-tririfirxf^r^r^f^^ 

I <■ CCCCCCCCCCCCCC-eCCOCCLCCCCwCCCCCCCCClCCCCCCCCCCC 

IU + + + + + + ♦ + + + + + + + + + + + + + «► + + + + + + + + + + ++ + + ++ + + + + + + + + 

I »- ircact>rircccN£a:ir^TK:ccr-ccrrr-r^Fr^r:cc»£ccra:^cv^ccrr-<rcciri s ‘Cir^r-c<vc\.irrr 

i k j fc vrr-x^c»o^-cc^-r-K!iru‘vr:c:cccr^X4ffr^-r-^cLr^iraK:r-vrvC*-cf ,> -crrrccor = vCcv 

i i/: vCrrLrr^s£f^a:^fCv-^cv.‘frf^ir. ^ovc^ca.^-r*'. r-*Mr-ccfCr-*H^v£c\r-r-cv.^-*^roif,tvf^sCfvr-\£ 

i h j. ^ ir c\jxO'X^x)irr^x.r»-t^tcv-- j ^rv,^^'v,^.^- — — Hri\nH 

I h 

I <x _J 
I K <r 
I If > 

I 

I 

I 

I lflfCHC‘£vO^lflflflf\C\CvCXvDvD\r;JJ'^d-^^Xf^^^^vrXvD\ 0 d-^^tfvfvfvCvC'f'f 

l 2 ccccczccccccc-ccc-c=rcrcccccccccccccccccccccc 

Lfl C C ♦ 4 - + I + + + + + + + + + + + + + + + -f + + + + -f + + ++ + + + + + + + + + + + + + + + + 

y i a*- if. cccr — cv^a: — ffLfs<f^hvCif^ — "c u ^ r** »c crKivf^^cv — cccxccr-rr — jrrrccarct — a c\ 

z i <thNC\£f^^o , ^a(vo i ffff'vC»rcccir. trecffCH^Hffd-cifcar.vfi-'vCcccK'. d-ifif^crf^a. w 
S I D<HOH IlO Hirn C Oi o.HwCr-HCMO'^Hr*- 

*- I Z hn ^ .H — tf CNJ —< <\J ~ -I — «->* fT- (V O' rC CV :t fT <V :t If C\ «0 K. CV (V C\. KV — K) fT 3 CV, CV CV lO <r- K rO :± fT- fO 



< I KLJ 
5 I If C 

I 

I 

I 

I 

I 

i x f : ^ c »c r £ f f f if, x v: x r >c x : i ^ i i x x -f x j t i i x -: £ x i ^ £ x ■£ r ■: x 

I c-ecceoccce-ececccc.-cs&cc = =.ccc.e.-©cccce&cc=cGcc©cc 

I + + + I + + + + + + + + + + ++ + + + + -f + + +*t- + 4-f + + + + + + + ++ + + + + + + + + 

I I h- if ^ ^ t C i {V c (O vf A’ N (V >£ T I'*' ? f>. c N f r* — i 1 rf ^ (V ? w T «-< c .V J ff. f 

i uJ , fia'roor*-'f>-HAjrr;a)^^X5Loz'c\j^j'©©v0^f>r^x-r^.Mcvrrr;c’3o = >cwf)c:X3ij'ro k 'va'fO 
I zz)X)XLr)xaDj-oj«\jcffOor , ^f s ~vOXJ-r-XLr.r-*xa:r-c(v^'d->wC? s -a'^zfAJpoc\jirxr s -r s ~j-o^cro s x 
I Lrr-rofo^in^irirh-frrod-^^ir. itn ^ ^ 

I It 

I 5 > 

1 

I 


H w H (\| ro f- A f H (V K) —• OJ K> n <Vi K, ^ CV rC .- CV K; — (V K) « Mf ^ (V f — i C\. K) CV 

L. I *.<1<«XXXX «C < < X X X <X <1 < « <X <. =S < <1 X X 

51 I K H” H H H H H H H H 

<x l U (jH^^^^^^^^^^^(vcvcvjaa<vwacvc\ ; "^^w^H^^r)K'fOrrK',fOirrfd- =t 

X I < 4 S y 5 5 > s -yn /1 / 4 r A S 5 S 5. 2 (/; /. »; H- s. >, k JS If V, O- J 5 £ if {/! I/ 1 . < r S i S i ?. 


' i ^Airiiffhco <T^<\rr^irxr-cri7c^<v»*~^fir«xr*-cr<rc^-cvp* r j‘tr4‘^cirrc — rvyrd-x 

Cl CvrvjfMCJCJK^fOfOrrrOfT'OfOrrPO^'^-d-d-cJ-d- 

£T I 
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TABLE 11. (Concluded) 



x, ir. x tr.vCf^f^r^irirtr^xCf'-r^r^r^r^r-iriririfvCvCNCNtriririnirvCxc^irtruiirx 
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC'CCCCCC 
+ + + + + + + + + ♦ + + + + + + + +♦ + + + + + + + + + ♦♦♦ + + + + + + + + + + 
crirrrcvcf^ircc.a»ro:*-icircr^cc — i^ircvvccC'^ccir. ff^Kccvcrf^^NCzri^vCK: 
cco , vccccvK:o'<vO'e-^*rc!y'3'r-0'irir^*r-c-'vCvcc\fvh-f*"itwc 
h c 10 ir. ^ (vmt. - jxcroa vC iocr a cvr^irr-c%cir>OK)oc*Hd-ccch'CNjO(c*-<\CO'c 
<\]^^^^-'X)^wcv^*H(\)cC\x:^*-*-'"^^-^r-icvj-rvC(T'vC^HK:f\iro<vjc3-^s0cc^«vjcv»H(\/ 


CCCCZCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC'CCC 
+ + + + + + + + + -f + + + + + + + + + + + + + + «. + + 4 . + +4 + + + + + + + + + + 

^^»rCT%rcc — — r^o*-ccK:<\*-^^r-(Tcr*c<\ir — ^acairo‘ccevtr^r^cj'sCscc<r 
i^K:ircCifw0‘r«‘N^r-K<rf^Ki^cirotr<rcr(rif<NJsc^r^<vj<cirir!^r-0'NCCK5(r^'>£> 

^ oj <v cv. »r n ir to c* *-< — x * 


c * * x * x x x x a * * * x x x -x x vr x d> * * * x x -X x * ^ ix * in x x x ^ sn 

ccccccccocccoococoooooccccooccoococcooooc 
♦ + + + + + + + + + + + + + + +■►+ + + + + + + ♦ + + + ++ +♦ + ♦ + + + + + + + 
r- e :t tv x *- \C c c « x h a ^irtr-«.-''Xif;*-occj'fXX{ur^«rc>wK)3‘Xftjexaxo‘f^<v 
Ki^^r^inc«CNja'iOui^cr^co^^r*3K)®^^30ccfC^-rocofvifOcoina'^incc^oco 
a: ro^od^^O'r^r^trctcxf'-incticti.HCEXvoor* cNcocjcoa'Hcc^oNxoHXioo 
XfOd-ir, dH^^x^irx(\;^^d-x^K)irjirx(\j^iOHro(MnHi/)H<\iH(\jxxccinx 


K> «-< (V tr- ^ CM fO <V »0 ^C\JlC^CVlCw<V»0 HWICHWK) ^ <\j »o «-» oj ro 

I <«« < < < I £ I CCXCC I I IE 

c\;c\j<\icvi<\;<vj<\Jc\i!OfOK}rOfOfO»OfOiororr> oooooooo 

s tr. .r tr. < *- ssis s. str tn tr ^ rssmww<*-xisw»AUKh 


XNtfCc^tvtr^irxf-xac^ftfrdifxhffac^ftjrsirxrfftrc^Mr^trx 
^ ^ a-^trtr>tXixtr>inir;n*./iirxxxxxxvC\XXxr*-f*r^Nr-r*Kr*^N«cc*®cococo 
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TABLE 12. STATISTICAL VALUES - OUTPUT SET B, ENGINE 2 _ START POSITION 


_J — < «-H HWWWOOCHftlWftJOCOrtWftlWHO 

« cc cccccccccccccccccccc 

O 4 4 + 44144444414 44144414 

ct cr e c* m x rr cccffifHCMrirKiyCvCFCtfac 

»- mcccii-ircMtrcrxx^accrxirx^xxmr-© 
cr ir^cKivCciatnA^fo^aiocwc^^^Kiir 

K »«*••«»••••»•«♦ ••« •««•«« 

< J I I | I I t I I I I I I I I I I I I I I I I 
t- <. 

IT. > 


COK)(\IWWHCCOlOHC\|HWOCOi(\|W^ <MOCO 
cccccccccccccccccocccc cccc 
44 1 444 I 4+.4 1+44 14 + 4 I 44. ♦ • + + 4 

cec^cr-cacM— 'XM— '<v^M»ciirxc % ^Mcc a- <r x x 
;x^ 3 rcv’ircMrr>xxx.-*cMMcvc:<r a: r- x in x m cv cr cc r- 

0 'C(\.CirK)HhC 3 ’ClTHirH<tr-CrtN(Vav£CCHC 
XXCMCM—‘-H-H»rX. 3 -rOffi~HMf r >~«.d-CM^^^C 0 h H ^ (\J 


t I I I I I I I I I I I I I I I I I I I I I I I I I 


\r 

2 : 

s 


cc — — —>rr — — — —— — CV — 
3 c o ccoccccccccc 
C C 4 4 + 441 I I 14441 I 

(THO'CClOCClOXi’d'ACC t Ch 

ch-\i:ccJ^^ o x rh - 3 c : ir 
ccj^fvc: xrrxxrxxMCT cr — — 
z*-irirc«'-'rrirccxxx— >cv tv a> 



K Ll 
LO Q 


— — — CV — CC :r — — — cv — c 

cccc CCCCCOCCCCOCO 
I I 14 4 4 1 I 4 4 I 4 4 4 I 1 + 

•3- — xr-xx:xcr— (j-j'cr-xxroro 

r- ^ a x c MifccircftircKr- 

(V — fTCC— : 3 CV — CCCV X.XCC 
<t \C r, »C rrK-, ^ N -H h ^ KHV r. CC 


— d- — — — — — '3 - — — — <y — — —< 

OOOOOOOOCOOOCD COCC 
-11+441 II I 14 4 4 I I II 
xr-crcccv^ird-cvr-ccxcr hcccn 
co (Vffifirr<v^ircrr <v r ^ r. ^ 
a — Mctd-NcrcrrorrMMx rr a a r. 

CtO'w(\|HO'^(£tfHl< 5 ftIr' WK 1 CS 


OO ^HHHlOHOCHHHlOrt 

cc ccccceccecce 

4 4 + 4411444++11 

x cv e crxirr-cv^c— xr-rr;r 

itcccK:cicc>C(vcMj'Hft!C 

rrcr*cifao :scccirr-r-irc. 

L. <x • • * • • • 

2 > I I III I I I I I I I 


o o x ^ iK> — o c ^ — i— * «h m c o 
ccccccccccccccccc 
+ 4* 14 4 4 1 I 4 4 I 4 4 4 t 4 4 

f^irccc^c^ircvcuirr^ccvffr^ir 
^N^CHCifina'cvuccmcffoc 
ftccK. ccif. 3ao'ca3c<rc«i 

— IT IT. 3* ^ ^ O'HHr^vC^^ 


II III III til It 


lO H O H 

ccccccccccccc cccc 

4 14 4 4 1 I 4 t 14 4 4 I 14 1 

(tcfftt^airactovccd' r. c ir cr. 
Htfcccc^Mrft!<cN^(re cm x cm x 
^•rir — K'lTst — cr — ir^r^ m ;r — c 
— i^cvprcv — x—txcM^rrcv ro x x 


i ill ill ill ill 


_I — — — 1 CV CV CV C C C — ♦ CVCVCVCCCCCVCVCVCVCCCCCVCVCVCVCeeC — CV — (V C C C C CV CV (V c c c 

c ccccccccccccccccccccccccccccccecccccccccccccc cccc 

{J 44 + 444444+4444444444 + 444444 + 4444444 + 444444444 4 4 4 4 

*- rr, oirf^irr-vc — c\ r- c cr cv — c\ x x <r — < cc ^ x ~ freer- — cc x c x cm — — r- cr m cr x rr xaiITMT cv — x it 

t- xrtxcv — -'Cr^fxxcxxacvcrxcvcccrcvxaccrcvccvc-r'-— xcr^ccxccv —icvcrcv— 'Xcv <T cv cv cr 
x u*j^C'rc''aCvca.« i r a < cc m x r- m m — .axe— icrccrrctMcr:*;* — i io a: s a c t ? n ^ n h a* cc 

— u_ — — CM cf — — cvcv-’T^X — — — m.'Ornj' — — — cvmxxx — — — (vrTzr^XX — M — — .4 — — X (\j — cf — 



1 

f- 

3 

• 

• 

« 

• 

• 

• 
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• • 
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<r 
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e — 
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•-H rH 

H 


«-< 

«H 


— « 



1 


Z 

c 

” 

c 


c 

c 

c 

c 


c 

c c 

c 

c 

c 

C 

c 

c 

c 

c 

c 

c 

c 

c 

3 


c 


c 

c 

c c 

c 

z c 

c: 

c 

- 

c c 

c 

c 

c 

C 

c 

C 

c c c 

X 

1 

c 

c 

4- 

4- 

4- 

4 

4 

4 

4 

1 

i 

1 

4 4 

4 

4 

1 

i 

1 

1 

4 

4 

4 

4 

1 

4 

4 

4 

4 

4 

4 

4 

1 1 

I 

4 4 

4 

4 

4 

1 1 

I 

1 

4 

4 

4 

4 

1 1 1 

z- 

1 

Ct 

*- 

a 

X 

Jt 

CT 

o 

X 

X 

c 

CV 

c 
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APPENDIX 


METHODOLOGY 


The development of the equations used to define the response- of the 
SSME for the side load excitation is presented herein. The formulation is 
structured to facilitate the application of the Fortran Matrix Analysis (FORMA) 
system for solution of the problem on the digital computer. FORMA is a library 
of subroutines coded in Fortran IV for the efficient solution of structural 
dynamic systems on an electronic digital computer. The FORMA subroutines 
are in the form of building blocks which can be put together to solve a large 
variety of structural dynamic problems. 

For convenience the methodology development is divided into two sec- 
tions: ( l) transient response analysis and ( 2) transient loads analysis. The 
equations of motion and the coordinates used to define the system state are 
explained in the transient loads analysis section. The transformation equations 
applied to obtain loads ( and displacements and accelerations) from the system 
state coordinates are explained in the transient loads analysis section. 


TRANSIENT RESPONSE ANALYSIS 


The transient response analysis involves describing the interaction of 
the structural model and the actuator models for an externally applied forcing 
function. The structural model is comprised of both the engine structure and 
the orbiter thrust structure. The actuators are simulated as a feedback control 
system, controlling the motion of the engine. The equations of motion for the 
transient response are obtained by: 

1. Developing the equations of motion for the structural model; replac- 
ing the actuators with externally applied forces. This matrix equation will be 
referred to as the structural equations of motion. 

2. Expressing the block diagram model of the actuators into a matrix 
equation. This will be referred to as the actuators equation of motion. 

3. Coupling the structure and actuators by applying compatibility 
relationships to obtain the equations of motion for the coupled system’s transient 
response. 
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Structural Equations of Motion 


The equations of motion for the engine without actuators are developed 
first. The actuators are replaced with external forces {F (t)}, (F, (t)}, (F (t)}, 

cl D C 

and{F ( j(t)} applied at the actuators attachment points a, b, c, and d. The side 
load forcing function (F e (t)j is applied at the centerline of the bell exit plane, 
point e. These forces are shown in Figure A-l where 


{F a (t» 


F (t) 

X 

F y (t) 

F z (t) 


3x1 vector of components of force acting at point 
a (same notation for other points). 


X 



Figure A-l. Applied forces for engine without actuators. 
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The discrete equations of motion of the engine without actuators included 
is written as 


[ M]{*h} + [D]{h> + [K]{h> 



F b («) 

F d (t) 

F a (t) 

F o (t) 


( A- 1 ) 


where 


i means a properly positioned matrix and 


{h} = discrete displacement, 

[M] = discrete mass matrix, 

[D] = discrete dumping matrix, 

[K] = discrete stiffness matrix. 

It is extremely convenient to transform the discrete equations into a 
more practical form. The mode shapes of the dangling engine ( engine without 
actuators) is available from the modal analysis. With these mode shapes one 
can write the transformation equation 


{h} = [<p]{ q} , (A-2) 

where 

[ <p ] = modal matrix, each column being a mode shape, 

{ g} = vector of modal coordinates ; size ( n x l) , 
n = number of modes retained in the transient analysis. 

The mode shapes are orthogonal such that 


[ 0 J T [ M] [0 ] = [Meq] , 


(A-3) 
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where ^MeqJ = diagonal matrix of equivalent masses. Additionally the modal 
vectors are energy normalized such that 


[Meq]=[l], (A-4) 

where [0 is the identity matrix. 

T 

Substituting expression (A-2) into(A-l) and premultiplying by [ 0 ] , 
one obtains 


[0 ] T (M1 [0 l(q > + [0 ] T [Dj [0 ]{q} + f0 ] T [K] [0 ]{q} = (0j T 


F e (0 


+ [01 T { 


F b (t) 

F d (t) 

F a (t) 

F c ( l ) 


(A-5) 


In view of expressions (A-3) and (A-4), the triple matrix product [0 ] [K] [0] 
reduces as follows: 


[0] T [K J [0 J = 



(A-6) 


where 



is the diagonal matrix of modal frequencies squared. 


One assumes that the damping matrix is such that 


[0] T [D][0] = 



(A- 7) 
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where 



diagonal matrix for model damping, 


= modal damping coefficient for mode i. 

With expressions ( A-3), (A-4), (A-6), and (A-7), one can write 
equation ( A-5) as 


+ 2 5 1 ^ 1 ]{q) + £ a. 2 j{q} = t0 1 J T {F e (t)> + ['? 2 ] T { 


F b (0 

F d (t) 

F a (t) 

F c (t) 


(A-8) 


where 

[0jJ = modal matrix whose rows correspond to the X, Y, Z displace- 
ments at point e; size (3 x n), 

[0 2 J = modal matrix whose rows correspond to the X, Y, Z displace- 
ments at points b, d, a, and c, respectively; size ( 12 x n). 

Expression (A-8) is the equation of motion in modal coordinates for the 
structural model with the actuators replaced by the applied forces; {F (t)}, 

cL 

(F, (t)}, (F (t)}, and {F ,(t)}. These forces are unknowns at this point in the 

DC Q 

analytical development and must be eliminated. 

Actuators Equations of Motion 

The block diagram representation of an actuator is shown on Figure 5. 
This block diagram represents the actuator grounded at one end with an external 
force applied at the free piston end. For the analysis here, the actuator con- 
necting points a and b are referred to as actuator no. 1; actuator no. 2 
connects points c and d ( see Figure A- 1). Let points a and c correspond 
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to the grounded ends of the actuators. It will be necessary to remove these 
’’ground" constraints since points a and c are points on the flexible thrust 
structure. 

Each block in the diagram contains the Laplace representation of a 
transfer function of the form 


= ^ E 
out D( s) in 


where D( s) is of quadratic or lower order and the order of N( s) cannot exceed 
the order of D( s) . This expression is written in matrix form as 


[D( s) -N( s) ] 


out 


E. 

in ; 


{ 0 } 


The several transfer functions involved in the block diagram can be 
stacked in one matrix expression of the form: 


[ H ] { E} = {0} 


(A-9) 


The blocks are connected by applying the compatibility relations which exist 
between various input-output coordinates in { E} as defined by the block diagram. 
These compatibility relations are expressed in a transformation equation as 


{E} = [CT]{6} 


(A- 10) 


66 


li 11 U 11 II II u u 


II II li u 11 


11 11 M 


HULL 



where 


[CT] = coordinate transformation matrix, 

{6 } = coordinates chosen to represent the system block diagram 

equations; size ( m x l) , 

m = number of block diagram coordinates. 

Substituting expression (A-10) into (A-9), one obtains a matrix equation 
for the block diagram as 


[G){6 } = {0} 


(A-U) 


where 


[G] = [ H ] [ CT ] 

G(I, J) = a.. S 2 + 13.. S + y. . 

ij ij iJ 


Then equation ( A-ll) can be expressed in a more convenient form as 


[a ] {6 > + [13 ]{6 } + [y]{6 } = {0} 


(A- 12) 


Now, a procedure by which a block diagram model can be expressed in 
a matrix equation which has the same form as the structural equations of 
motion ( A-8) has been outlined. This procedure was applied to the block 
diagram of the actuators. The equations for both actuators are stacked in one 
matrix equation of the same form as equation (A- 12) . The bypass load limiter 
loop is omitted for this analysis. The matrix expression obtained for the 
actuators equations of motion is as follows: 
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W 


•• 

{«} 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


f fd) 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


7(2) 

0 

0 

Mp(l) 

0 

0 

0 

0 

0 

0 

0 


Xpd) 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 


Q v d) 

0 

0 

0 

0 

A(l)T p (l) 

-K p (l)T p (l)K o (l) 

0 

0 

0 

0 


<v» 

0 

0 

0 

0 

0 

A(l) 

0 

0 

0 

0 

Xj(l) 

0 

0 

0 

0 

0 

0 

M p (2) 

0 

0 

0 


X p (2) 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 


Qv(2) 

0 

0 

0 

0 

0 

0 

0 

0 

A(2)T p (2) 

-K p (2) T p (2)K o ( 2 ) 


?dpf<» 

_° 

0 

0 

0 

0 

0 

0 

0 

0 

A(2) 


Xj(2) 


0 

0 

0 

0 

Vi>y i)K o (1 > 

0 

0 

0 

0 

0 


0 

0 

0 

2f v (l)w v (l) 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

A(l) 

1 

0 

0 

0 

0 


[ 0 ] 

0 

0 

-K o (D 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

K p (2)T p (2)K 0 (2) 

0 


0 

0 

0 

0 

0 

0 

0 

2f(2)w v (2) 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

A(2) 

1 


0 

0 

0 

0 

0 

0 

-K q (2) 

0 

0 

0 


{«} 


* 1 

f(l) 


f(2) 


Xpd) 


Oyd) 

Xj(l) 


X p (2) 


Q»(2) 

^dpf 

Xi ( 2) 



0 0 


-1 

0 

0 

0 


0 

0 

K 0 (l) 

H(1)K v (1)u; v 2 (1) 

0 

0 

0 

0 

0 

0 


0 

0 

0 

“ v 2 (l) 

0 

-1 

0 

0 

0 

0 


m 

0 0 
0 0 


0 

0 

0 

0 

0 

0 

K o (2) 

H(2)K v (2)u> v 2 (2) 

0 

0 


0 

0 

0 

0 

0 

0 

0 

w v 2 (2) 

0 

-1 


0 0 
0 0 


0 

0 

0 

0 

0 

0 


{ 6 } 

f(l) 

f(2) 

Xpd) 

Oyd) 

DUMMY 

DUMMY 

X p (2) 

Qy(2) 

DUMMY 

DUMMY 
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In the above expression the index ( l) indicates that the parameter is for 
actuator no. 1 and an index ( 2) indicates actuator no. 2. The block diagram 
parameters are described in Table 4. 


Coupled Actuators/ Structure Equations of Motion 

In the previous sections the equations of motion have been developed for 
( l) the structural model without actuators and (2) the actuators with one end 
grounded, the other end free. Now these systems will be coupled and the equa- 
tions of motion for the combined system will be obtained. In order to accomplish 
this, compatibility conditions between the two systems involving force relation- 
ships and displacement relationships must be enforced. 

One knows that, if one neglects the mass of the actuator, the forces 
applied at the attachment points are opposite in direction and equal in magnitude. 
Then for actuator no. 1, 

F = -F b - f( 1) U, , 


where Uj is a unit vector along the actuator axis directed from point b toward 
point a. Similarly, for actuator no. 2, 


F c = -F d = f(2) U 2 , 


where U 2 is a unit vector along the actuator axis directed from point d toward 
point c. 


One can express the above force compatibility conditions in a matrix 
form as 


{ 

{ 

{ 

{ 




[ TFf] {f} 


(A- 14) 
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where 


[TFf] 


-{Uj 

{0} 

{0} 

-{U 2 } 

{Uj 

{0} 

{0} 

(U 2 ) 


and 





Substituting equation (A- 14) into equation (A-8), one obtains the expres- 


sion 


[i]{« * [ 2^.J{q} + = [0l] T (F e (t)> + [0 2 ] T [TFf]{f> .(A-15) 


Since the actuator forces, {f}, are unknowns, one chooses to rewrite 
equation (A-15) with the actuator forces on the left-hand side to obtain the 
following expression for the structural equations of motion: 


[i](0] 


{q} 


h][« 




* + 

L ii 





[0] [0]_ 



pi J |p 2 l T fTFf]J 

|{q} 1 

[0] [0] 

|{f> 


[0J 

[ 0 ] 




(A- 16) 
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Now, the structural model equations of motion, (A- 16), and the actuator 
models equations of motion, (A-12), can be stacked in one matrix expression 
as follows: 


'[Olof 

{q} 

► + 

M 

[of 

, 

{q} 

> + 

’ Kh 1 ’ 

(q> 

^ [ 0] [a ]_ 

>{«} , 


[0] 

63]. 


,{5} 


[0] [y] _ 



[D,]{F e (t)} 


(A- 17) 


where 


T 

[<Pl] 

1 

[ o] 

] 

_ _ 

A 

m — 3 — ► 



and 


[D 2 ] = 


m 2 

*(m-2)* 

_[-[0 2 ] T lTFf]_ 

[0] _ 


I- 


The force compatibility conditions, (A-14), are contained in equation 
( A -17) ; however, there are displacement relationships yet to be included. The 
ground constraints for the actuators are relieved and deflection compatibility 
conditions are implemented by a transformation which will be applied to equation 
(A- 17). 

Consider actuator 1 only for the time being and let the vectors p and 

3. 

denote the displacements of ends a and b, respectively. These vectors 

are shown on Figure A-2. Then, take components of p and p, along the 

a b 

undeformed actuator as follows 
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a 


6 


(P a ' U t ) U t 


and 

? b = . U^U! . ( A- 18) 



Figure A-2. Actuator geometry. 


One assumes that the deflections will be small and relieve the ground 
constraint at point a by applying the expression 


6, = 6 + X (1) Uj 

b a p' ’ 1 


(A- 19) 


Substituting expressions (A- 18) into (A- 19) and solving for X ( l), one obtains 

P 
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(A-20) 


X p (i) = 5 ■ 0> b -P a ) ■ 


The displacements at the actuator ends are equal to the displacement 
of the structural model ( engine and backup structure). Using the mode shapes 
of the structural model to define the deflections of the structural model, one 
applies equation (A-2) and expresses the displacement compatibility conditions 
as: 


a 


-W> 


t^ a J {q) 


and 


- W 


U b ){q} 


( A-2 1) 


where 


[$ ] = modal matrix whose rows correspond to the X, Y, Z displace- 

£L 

ments at point a, 

[0, ] = modal matrix whose rows correspond to the X, Y, Z displace- 
b 

ments at point b. 

Substituting expressions (A-21) into (A-20), one obtains 

X (1) = [TXPlQ]{q) , (A-22) 

r 


where 


[ TXP1Q] = {UJ T ([? b l - [<P a ]) 
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In a similar manner, the following expression is obtained for actuator 2: 
X (2) = [TXP2Q]{q> , (A-23) 


where 


[TXP2Q] = {U 2 } T ([0 d J -itj) 


Using equations ( A-22) and (A-23), one can write the following trans- 
formation: 


{q> 



co ! 





0 



] 


r -n 

to} 

f(D 



0 

"1 

1 

1 

0 

0 

0 

0 

0 

0 

o' 


f( 1) 

f(2) 



0 

1 

1 

1 

0 

1 

0 

0 

0 

0 

0 

0 


f(2) 

x p (i) 



ITXPIQ] 

1 

1 

0 

0 

0 

0 

0 

0 

0 

0 


« V U) 

« v (» 



0 

1 

1 

0 

0 

1 

0 

0 

0 

0 

0 


Dummy 

< Dummy 

►s 


0 

J 

1 

1 

0 

0 

0 

1 

0 

0 

0 

0 


Dummy 

Dummy 



0 

1 

1 

1 

0 

0 

0 

0 

1 

0 

0 

0 


Q v (2) 

X p (2) 



J_TXP2QJ 

I 

1 

1 

0 

0 

0 

0 

0 

0 

0 

0 


Dummy 

Q v (2) 



0 

1 

1 

0 

0 

0 

0 

0 

1 

0 

0 


Dummy 

L J 

Dummy 



0 

1 

0 

0 

0 

0 

0 

0 

1 

0 



Dummy 

L J 



0 

1 

1 

0 

0 

0 

0 

0 

0 

0 

1 




or, more concisely, 
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(A-24) 


{q> 


w 


(q) 

(0. 


T 

Substituting equation (A-24) into equation (A- 17) and premultiplying by [tp ] , 
one obtains the equations of motion for the combined structural and actuator 
systems as 


[A*] 


(q) 

• • 

u> 


+ [B*] 


{q» 

{«} 


+ [C*] 


(q>. 

= [d*]{f (t)} 

{«}[ 


(A-25) 


where 


[a*] = 


[IJ [o] 

[o] fa] 




[B*] = [^] J 


[ 2 $ w i] [ 0] 

[o] m 


[4>\ 


[c*] = [ipY 


’H] tD 2 ] 

. [o] r r IJ 


[4>] 


9 


[D*] = I^J T [D t ] 
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It can be shown that the rows and columns of [A*], [B*], and [C*] 
which correspond to f( l) and f( 2) must necessarily be zero. Likewise the 
rows of [D*] corresponding to f( l) and f(2) must necessarily be zero. [The 
coupling procedure could have been accomplished by solving equation (A- 13) 
for f( l) , f( 2) and inserting into equation ( A- 15) , thus eliminating f( l) , f( 2)]. 
As an error check, the rows and columns corresponding to f ( l) and f( 2) 
are examined for nullness in the computer program. By collapsing the null 
rows and columns from the matrices, one obtains the equations of motion for 
the transient response as: 


[A] 


{q} 


N) 



■+ [B]i 

+ [Cl 

(p) 


|{P}J 

Up)] 


[Dj{F e (t)> 


where 


(A-26) 
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X(D 
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Q v (2) 


Q v (2) 


Q v (2) 

V 2) 
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Dummy 
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x ( 2) 

1 


Dummy 


and [A], [B], [C], [D] are obtained from [A*] , [B*], [C*], [D*] by remov- 
ing the null rows and columns corresponding to f( l) , f( 2) . 

The transient response equation (A-26) is solved numerically for 


(q) | 

|{q) 

(p}j 

y { 

UW I 


and 


(q) 

{p} 
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at incremental times using FORMA subroutine TR2. This subroutine uses a 
third-order Newmark- Chan- Beta numerical integration procedure. 5 The 
transient response solution is saved and used in load transformation equations 
to obtain the transient loads experienced by the system. 


TRANSIENT LOADS ANALYSIS 


The equations used to obtain internal loads ( shear, bending moment, 
torsion, etc.) are developed herein. A loads transformation matrix, [ LTM] , 
is available from the vibration analysis. This matrix transforms the total 
force acting at some time, t, to the internal loads, thus: 


{ L (t) } = (LTM]{F(t)} 


(A-27) 


where 

{ L(t)} = vector of internal loads, 

[ LTM] = force loads transformation matrix, 

(F(t)} = total force vector, (includes inertial and applied forces). 

We know that the total force acting on a structural system can be 
determined from the displacement as follows: 

{F(t)> = [K]{h(t)> . (A-28) 


5. Chan, S. P.; Cox, H. L.; and Benefield, W. A.: Transient Analysis of 
Forced Vibrations of Complex Structural-Mechanical Systems. J. Roy. Aeron. 
Soc., July 1962. 
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Substituting equation (A-l) into equation (A-28) yields 


(F(t)} = - [ M] { h (t) } -[D]{h(t)} + {|F e (t)} + 


< F b( l > j 

! F d (t) ! 

Im4 



(A-29) 


Now, neglect the dissipative forces (generally conservative assumption) and 
substitute equation ( A-2) into equation ( A-29) and obtain 


{ F(t )} = -[M] ]{q (t)> + {|F (t)|> + { 

i C | 


,F b (t) 

F d (t) 

Fjt) 

F c (0 


(A-30) 


Substituting equation (A-30) into equation (A-27) yields the following expres- 
sion for internal loads: 


(L(t) } = -[LTM] [M]l^]{q(t)} + 


[LTMl]{F e (t)} * [LTM2]j 


F b (t) 
F d (0 , 

F a<‘> 


(A-31) 


F c (t> 


where 
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[LTM1 ] = submatrix of [LTM] whose columns correspond to X, Y, Z 
forces at point e, 

[LTM2 ] = submatrix of [ LTM] whose columns correspond to X, Y, Z 
forces at points b, d, a, and c, respectively. 

Substituting equation (A- 14) into equation (A-31) yields 


{ L(t) } = -[ LTMMP]{q( t)} + [LTMl]{F (t)} + [LTM2] [TFf ]{f(t)} 

© 


(A-32) 


where 

[ LTMMP] = [LTM] [M] [0] 

Now referring to Figure 5, the actuators forces f ( l) , f( 2) can be written as 

f( 1) = M ( 1) X ( 1) - K ( 1) X.( 1) + K ( 1) X ( 1) 

P P o 1 ° P . (A-33) 

f( 2) = M p (2) X p (2) -K q (2) X.(2) + K q ( 2) X p (2) 

We know from equation ( A-26) that 

X.(l) = p(3) 

X.(2) = p(6) 

Then using equations (A-22) and (A-23), one can write expressions (A-33) in 
matrix form as follows: 
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f(l) 

(f(t)} = 

f (2) 

M ( 1) 0 

P 

0 M (2) 
P 

[txpiqJ 

[TXP2QJ 

(q) 


X 

-K 0 (l) 0 


'm 

1 -t- 

K 0 (l) 

0 

[txpiqJ 

T 

0 

- K o (2) _ 


«6) 


0 

K q( 2 ) 

|_TXP2Qj 


(q) 


. (A-34) 


A more convenient form of equation (A-34) is expressed as follows: 


{«t)> = la] 


{q>. 


w| 

{q}] 


'+ lb]\ 

' ! 

l j 

■ + [c] 

.{p}, 


!<p>] 

(P> 


(A-35) 


where 


[a] = 


M p (l)|_TXPlQj | 0 0 0 0 0 

M (2)|_TXP2QJ j 0 0 0 0 0 


lb] = 


0 j 0 0 -K (1) 0 0 0 0 

o v 


0 0 


0 0 -K 2 

o 


[c] 


K Q (l) LTXP1QJ [ 0 0 

K q ( 2) |TXP2Qj ! 0 0 


0 

0 


0 0 0 

0 0 o 


Substituting equation (A-35) into equation (A-32) yields the following expression: 
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(L(t)} = - [LTMMP]{q (t )} + [ LTM1 J{ ^(t)} + [ LTM2 ] [TFf] [ a] 

{q}| ({q}) 

+ [ LTM2] [ TFf] [ b] \ } + [ LTM2] [TFf] [ cK } . 

{P>) i(p>) 

Finally, one can express equation (A-36) in a form convenient for solution 
using the FORMA library as 

{q}| |{q) | ({q} 

{ L (t )} = [AA]{ } + [BBK }+ [CCK > + [DDj{F } , (A-37) 

(P}J |(P> ) |{P>) 

where 

[AA] = ([-[LTMMP]; [0]]+ [ LTM2] [ TFf] [a ] ) 

[BB] = [ LTM2] [TFf] [b] 


[CC] = [LTM2] [TFf] [c] 


[DD] = [ LTM1] 

The solution of the transient response equation (A-26) applied in 
equation (A-37) yields the transient-internal loads. Acceleration and dis- 
placement responses can be obtained directly by solving equation (A- 2). That 
is, displacements are computed as 

(h(t)} = [0] { q(t )} 

and accelerations are computed as 

{h(t)> = I0j{q(t)} 
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{ q } I 

{p }/ 
(A-36) 
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